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^ ■ Abstract 

. We consider a neutrino source at rest and discuss a condition for the existence 

, of neutrino oscillations which derives from the finite lifetime ts of the neu- 

trino source particle. This condition is present if the neutrino source is a free 
Q^' particle such that its wave function is non-stationary. For a Gaussian wave 

p) • function and with some simplifying assumptions, we study the modification of 

1^ ■ the usual oscillation probability stemming from ts- In the present accelerator 

experiments the effect of ts can be neglected. We discuss some experimen- 
tal situations where the source lifetime becomes relevant in the oscillation 
^ ' formula. 
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I. INTRODUCTION 



Neutrino physics is one of the most active fields in particle physics nowadays. Apart 
from the impressive results of the underground experiments concerning atmospheric and 
solar neutrino measurements, also reactor and accelerator physics have a large share in the 
evolution of this field and the understanding of the neutrino mass spectrum and mixing 
matrix (for recent reviews see, e.g., Refs. ll|J^). Since the phenomenon of neutrino oscilla- 
tions is now close to being an established physical reality and constitutes the most tangible 
window for physics beyond the Standard Model, it is very important to know if there are 
any limitations to the validity of the usual formula for neutrino survival and transition prob- 
abilities with which the experimental results are evaluated. Such questions have been 
discussed extensively in the literature in the context of the wave packet and field-theoretical 
approaches (see, e.g., Ref. ||^ for a list of references). 

In this paper we use the field-theoretical approach in the spirit of Ref. ||^, which has 
proved to be a general and unambiguous method to analyse neutrino oscillations, and con- 
centrate on the effect of the finite lifetime of the neutrino source particle. We study the 
problem under the following three assumptions: 

i. the neutrino source is at rest, 

ii. the source particle is not in a bound state and, therefore, is described by a non- 
stationary wave function, 

iii. the wave function of the detector particle is stationary. 

The first assumption is of technical nature and allows us to use the methods and results 
of Refs. and can probably be overcome f^, however, the second one is essential to the 
discussion in this paper. The third assumption simply means that the particle with which 
the neutrino reacts in the detector is in a bound state and the wave function of the detector 
particle does not spread with time [0. Note that we also neglect for both the production 
and the detection process any possible interaction with the background. 

The finite lifetime = l/F of the neutrino source particle has an impact on the neutrino 
oscillation probability in the following two ways: 

(a) a suppression of the probability amplitude as a function of F 0, an effect independent 
of L, the distance between neutrino source and detection, 

(b) a damping of neutrino oscillations through the coherence length L™^ caused by the 
finite lifetime and given by the factor exp(— L/Lp*^) |0 in the oscillation probability. 

The main topic of this paper is the derivation of a correction to the usual neutrino oscillation 
probability according to point (a) and a numerical study of this correction. It has been 
shown in Ref. |^ that the effect (a) of F is only present provided assumption ii) is valid. 
This seems to be the case in the important experiments of the LSND M and KARMEN fTT 
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Collaborations, which study z/^ z/g oscillations with /i"*" decay at rest as z/^ source.[] In the 
following we will use the source and detection reactions of LSND and KARMEN as a model 
for our field-theoretical treatment of neutrino oscillations. 

The paper is organized as follows. In Section II we derive the transition probability for 
~^ oscillations and show how the finite source lifetime changes the standard formula of 
neutrino oscillations and gives rise to a suppression factor for the interference terms in the 
oscillation formula. We will see that the corresponding condition for the existence of neutrino 
oscillations, i.e., the condition that the suppression does not take place, depends only on 
parameters supplied, at least in principle, by the experimental set-up and the neutrino mass 
squared difference Am^. This is natural within the field-theoretical treatment which enables 
the study of the dependence on those quantities which are really observed or manipulated 
in oscillation experiments. In Section III we apply the results of Section II to the case 
of a Gaussian /i"*" wave function and make a numerical study of the 2-flavour transition 
probability to investigate quantitatively the influence of the correction factor. In Section IV 
we present a discussion of our findings and draw the conclusions. 

II. THE CROSS SECTION 

As mentioned in the introduction, for definiteness and also for comparison with the 
notation in Ref. we discuss the process 

/i ^ + z/e + z/^ Ve + P n + , (2.1) 

which is investigated in the LSND and KARMEN experiments. As shown in Ref. 0, with 
the help of the Weisskopf-Wigner approximation one can write the amplitude of the process 
in the limit t ^ cxd as 



(2.2) 

where T is the time-ordering symbol and V is the total decay width of the muon. 'K^i^t 
Ti.^ jj^^. are the relevant Hamiltonian densities (in the interaction picture) for the production 
and the detection of the neutrinos, describing muon decay and proton to neutron transition 
(^.1[), respectively. The indices S and D denote source and detection, respectively. The 
muon /i"*" and the proton p are localized at the coordinates xs and xn, respectively. The 
proton state is stationary whereas the decaying muon at rest is described by a free wave 
packet with an average momentum equal to zero according to the assumptions i)-iii). Hence 
the spinors of the initial particles are written as 



^In the LSND experiment also Ue oscillations are investigated where the muon neutrinos 

originate from tt^ decay in flight. According to assumption i) we do not discuss this case here. 
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/oo p poo P 

dh J d^x^ J dt2 J d^X2nl,^,ixi)e-'2^''n-^,^^{x2) 



^,{x)=ijp{x-XD)e-'^-' with ^^{x-xd) = j,^^ j d'pi,p{p)e'^-^^-^^^ 

and ^i)p{p) = ^p{p)Up{p) (2.3) 

for the proton and 

^^^^M(Ple-^^'''"''^^'^*^><e^''''' with Mp) = i^',{p)v,ip) (2.4) 

for the muon with E^{p) = ^Jvn?^ + p^. In the 4-spinors Up{p) and f^(p) we have left out 
the polarizations of the proton and muon, respectively, because they are irrelevant in the 
further discussion. ^ 

The function is peaked at y = and the wave packet '^^{p) in momentum space 

is peaked around the average momentum (j?) = 0. The final particles are described by 
plane waves. After carrying out all integrations, the leading term of the amplitude in the 
asymptotic limit L ^ oo can be written as 

A^ = Y, U^jU:/''^''A^Af , (2.5) 
j 

where Aj and Aj^ denote the amplitudes for production and detection ( |2.1| ) of a neutrino 
with mass rrij, respectively. The parts of these amplitudes which are important for the 
further discussion are given by 

and 

Af = ---M-(l/+P2). (2.7) 

For the full expressions see Ref. 0. The kinematical quantities occurring in Eqs. ( p.6|) and 
( p.7[ ) are defined by 



Ed = E'^ + Kd - Ep , 

Es, = E^{q/+p^)-El-E'^s (2.8) 
Pi=p'u+p'es, P2=Pn+p'eD and L = \xd-xs\. (2.9) 



and 



Note that the first two formulas in Eq. ( |2.8[ ) follow from assumption iii) in the introduction 
@J^. Eq. ( p.6|) shows the dependence of A°^ on V (see point (a) in the introduction). The 
structure ( |2.5| ) arises from the fact that in the limit L — oo the neutrinos with mass m^, 
described by an inner line of the Feynman diagram derived from A ( p.2|) , are on mass shell. 
Consequently, the amplitude for neutrino production and its subsequent detection factorizes 
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for each j into a product of production amplitude and detection amplitude 0. Looking at 
Eqs. (|2.5| ), ( p.6|) and (2^) it is evident that oscillations involving rrij — ml can only take 



place if the conditions 

\Qj-Qk\^(^s and \qj - Qkl < (^d (2.10) 

and 

\Esj-Esk\<lr (2.11) 

hold, where as and cr^ are the widths of ip^^ and ipp, respectively. We call conditions ( |2.10| ) 
amplitude coherence conditions (ACC) and equation (|2.11| ) the source wave packet - finite 
lifetime condition (SFC) If they are not fulfilled then some terms of the amplitude and 
consequently the corresponding interference terms in the cross section (oscillation probabil- 
ity) are suppressed. 

The cross section a is obtained by taking the absolute square of the amplitude (p. 51), 
integrating over the final state momenta pi, pes, p 'n, pLd and averaging over the spins of the 
initial muon and proton. Hence 

a = J dP \A°°W (2.12) 

where dP denotes the integration over the momenta of the particles at the external legs and 
is given by 

,p _ d% d^p'^s d^p'^ d'^p'^jj 

and the integration is done over some volume of the phase space. 

In general, the integrations in ( |2.12| ) cannot be performed without knowledge of the 
source and detector wave functions. However, as noticed in Ref. 0, for 

r < os^D (2.14) 

the factors 

{[i{Es, - Eo) + r/2) X {-t{Esk - En) + r/2)}-' (2.15) 

in the cross section are strongly peaked with respect to Eo and we interpret the conditions 
(|1|) that the rest of the cross section is flat with respect to Ed if varied over intervals 
several orders of magnitude larger than T. This is a reasonable assumption because in the 
LSND and KARMEN experiments the stopped muons have momenta of the order 0.01 MeV 



TT,T2|, and thus as will be in the same range. Assuming atomic dimensions of the spread 
of the detector particle wave function, we find ao ~ 10~^ MeV. In any case, even if our 
guesses for as and an are wrong by several orders of magnitude, these widths can never 
be so small like the decay constant of the muon F ~ 3 x 10~^^ MeV. Therefore, we adopt 
the procedure that integrating over momenta of the final state of the detector leads to an 
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integration in the variable Ed and in view of the expressions ( p.l5| ) it suffices to pick out an 
integration interval of length Ai?£) which fulffils F ^ -C (Js^d- Considering neutrino 

oscillations with respect to — m\ = Am'^i^, then this integration interval should contain 
Esj and Esk- Since the condition ( p.ll|) is necessary for these oscillations to happen at all, 



we define a mean value 

En = {Esj) (2.16) 

for the neutrino mass eigenfields participating in the oscillations and AE^ only has to be a 
few orders of magnitude larger than T. In addition, the expressions ( p.l5| ) suggest to close 
the integration over the interval of length AEd along the real axis by a half circle in the 
complex plane and to apply Cauchy's theorem [j^. 

In this way we have shown that the cross section contains the damping factor 

leading to a suppression of the corresponding interference term in the oscillation probability 
when L becomes greater than the coherence length Lf^ = AE'jy/ Am^-fT . As already dis- 
cussed in Ref. [^, this coherence length is of the order of 100 light years for typical input 
relevant for the LSND and KARMEN experiments and its effect is certainly negligible for 
all experiments with terrestrial neutrinos. In the following computations we will neglect the 
coherence length arising from the finite lifetime of the neutrino source. 

Then, the only relevant integrand with respect to Ej:, is given by the factors ( |2.15| ) and 
we get 

1 1 _27r 1 

oo t{Es, - Ed) + |r -t{Esk -En) + \T 11 + t^{Esj - Esk) ' ^ ' ' 

On the right-hand side of this equation we can use assumption i) (the muon at rest) to 
compute to a very good approxmimation the energy difference 

Am% ^ _ ^ 

Esj - Esk ^ ^ l-{p{ + EdI). (2.19) 

As the next step, we want to perform the integration over d^p'^. Note that the momentum 
pi of the z/g from the source reaction (see Eq. ( p.l| )) cannot be measured. Again we make 
a simplifying approximation based on assumption i), namely that the muon wave function 
can be represented by 

Mp, + q,l) - ^'M + Ed) v^{Q) . (2.20) 



Then the integration over d^p'^ concerns only the functions {Esj — Esk) in (|2.18| ) and iIj'^. 
Therefore, we define the quantities 



g,k = I d\ W^m- ; , (2.21) 

1 - ipjkl-u/as 
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where we have defined u=pi + EJ and 



p.k = • 2.22 



Note that in view of Eqs. (|2.8| ) and ( 2.16| ) we identify the neutrino energy which can in 



principle be measured by the detector as 

E, = Ed. (2.23) 



The quantities ( |2.21|) have the properties that Qjk = glj and Qjj = 1. Since we consider 



ultrarelativistic neutrinos, we can take the hmit rrij — > Vj in all terms of the cross section 
( 2.12| ) except the oscillation phases and gjk- Then the probability for z/^ — > z/g oscillations is 



given by 

p- - = V If/ -Plf/ -P + 



g U:^U,,U,,U;, gj, exp (^-z^^ j j . (2.24) 

In contrast to the formula derived within the standard treatment of neutrino oscillations, 
the quantities gjk appear as correction factors in Eq. ( |2.24| ). These non-standard factors gjk 



represent the quantitative effect of the SFC condition ( 2.11|) discussed in Ref. [0. They arise 



from assumption ii) (see introduction) that the muon is described by a free wave function, 
which does not have a sharp energy, and would disappear for a muon in a stationary (bound) 
state (see Ref. 0). 



III. THE CASE OF A GAUSSIAN MUON WAVE PACKET 

In order to get a feeling for the effect of the (^-factors ( |2.21| ) on the oscillation probabilities 
we confine ourselves now to the case of a Gaussian muon wavefunction 

= iV^cTs)-'/' exp (-^) . (3.1) 



's 



Then we get 



gjk = / du e-^ I'^s , , (3.2) 



where u = l-u. Since the exponential factor in Eq. ( p.2|) is an even function with respect 
to the variable u the ^f-factors become real and defining y = u/as one obtains g^k = g{pjk) 
with 



g{p)^— dy ^ =^exp -U>J- , (3.3) 

^ J-oc l + p^y^ p VP / VP 
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where $c is the complementary error function (see, e.g., Ref. [|T5|). Note that g{p) ^ 1— 



for p -C 1 and g{p) a/tt/p for p — >• oo. Under the assumption of oscillations between two 
neutrino flavours the transition probability becomes 

■ sin^ 1Q\\- gi2 cos , (3.4) 



where 9 is the 2-flavour mixing angle. In comparison with the standard formula the factor 
5^12 in ( |3.4| ) lowers for a given transition probability the upper bounds on Am^ and sin^ 26. 
This can be seen explicitly by considering the low Am^-region of the parameter space where 
the cosine in expression ( |3.4| ) can be expanded in Aim?. If now p ^ 1 holds for an average 
Am?, in the low Am^-region we get p -C 1 and can use the corresponding expansion of g{p) 
given above. Then, in this region the isoprobability contour is determined by the relation 

IfiE^ P- - 

(ArnTsin^2g^ ^r;SCr. - (3-5) 



It is clear from Eq. (|3.5|) that a large source lifetime and a large source momentum spread 
lower the upper bounds on Aw? and sin^ 29 compared to the standard oscillation formula. 
For the source lifetime effect to be significant in a given experiment, the numerical value of 
<7sTs/^fi must be comparable to the source detector distance L. 

In the following we investigate quantitatively the effect of gu on the transition probability 
( W^ - In Fig- 1 we display the function g{p). One can see that at p = 0.5 the function g is 
around 0.9 which means that gi2 starts to deviate appreciably from 1 when p becomes greater 
than about 0.5. In the case of the LSND and KARMEN experiments, using F ~ 3 x 10~^^ 



MeV for the decay width of the muon and the typical numbers as — 0.01 MeV |[Tl| , p!2 
Am^ ~ 1 eV^ and E^y ~ 30 MeV, we obtain pi2 — 0.5 x 10~^. From this estimate we conclude 
that with the above input numbers the correction to the transition probability is negligible 
for the LSND and KARMEN experiments because l—gu ~ 10~^ (see remark after definition 
of g ( |3.3| )). Note, however, that pi2 — 0.5 x 10~^ is only two orders of magnitude away from 
having a 10% effect. This motivates us to have a look at the influence of pi2 {gu) on exclusion 
curves obtained, for instance, in the KARMEN experiment. Therefore, we use L = 17.7 m, 
average in Eq. ( p.4| ) over the energy spectrum of z/^ and as a Gedankenexperiment we vary 
the momentum spread as- In Fig. 2 we show the corresponding exclusion curves in the Am^ 
- sin^ 29 plane for the cases of cr^ = 0, 1 and 10 MeV corresponding to gi2 — 1, 0.9 and 0.3, 
respectively. We can see that for as = 10 MeV the exclusion curve has changed noticeably. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper we have studied the effects of the finite lifetime of the neutrino source 
particle. In particular, we have eleborated the condition for the existence of neutrino oscil- 
lations arising in a situation where the neutrino source particle is not described by a bound 
state and thus represented by a non-stationary wave function. Furthermore, in our analysis 
we have made the assumption that the neutrino source is at rest. In order to perform an 
analysis closely related to experiment we used the field-theoretical approach which allows 



8 



to work directly with the states of the particles responsible for neutrino production and 
detection and where the oscillating neutrinos are represented by an inner line in the com- 
plete Feynman diagram containing the source and detection processes. Thus the neutrino 
oscillation amplitude is determined simply through the neutrino production and detection 
interaction Hamiltonians. For simplicity, in the following we will concentrate on 2-flavour 
oscillations. 

The finite lifetime of the neutrino source leads to a finite coherence length |]T3[. In our 



case we have explicitly indicated the damping factor ( p.l7| ) from where the coherence length 



associated with F can be read off. Using numbers taken from the LSND and KARMEN 
experiments for the quantities appearing in we get a length of the order of 100 light 
years [0 which is absolutely negligible compared the the coherence length stemming from 
the fact that the neutrino energy Ed (see Eqs. ( pl8| ) and ( |2.16| )) determined by the energies 



of the final states in the detection process is not accurately known. Therefore, the effect (b) 
of F (see introduction) is irrelevant. 

Let us now compare the amplitude coherence condition (ACC) with the source wave 
packet - finite lifetime condition (SFC). The ACC can be reformulated as |l^,p|,|6| 



<J.S,ccD < ^L°^' , (4.1) 

where (Jxs,xD are the widths of the source and detector wave functions, respectively, in co- 
ordinate space and L"*^" = ATiE^/Am^ = 2.48 m MeV)(l eV'^/Am^) is the oscillation 
length. For the LSND and KARMEN experiments with E^ ~ 30 MeV and Am^ ~ 1 eV^ 
as found by the LSND collaboration, one gets L°^'^ ~ 75 m. This is to be compared with the 
widths crxs,xD which are certainly in the range of atomic distances or smaller. Therefore the 



conditions (4J-) are very well fulfilled. We want to stress once more that the SFC and the 



correction factor gi2 ( |2.21| ) corresponding to effect (a) of F (see introduction) derive from two 



ingredients: from the finite lifetime of the muon taken into account through the Weisskopf- 
Wigner approximation and from assumption ii) in the introduction that the wave function of 
the muon is non-stationary. The first ingredient leads to the factor l/{i{Esj — E^) + ^F) in 
the amplitude instead of the familiar delta function. The second one leads to the momentum 
dependence of Esj as given in Eq. (p.8|). If the neutrino source particle were in a bound state 
and, therefore, described by a stationary wave function, Es would be fixed and independent 
of J and momenta, and there would be no condition ( p.ll| ). Note that in our approximation 



the Hamiltonian which determines the wave function of the source particle does not include 
weak interactions. However, the decay of the source particle via the weak interactions - 
whether its wave function is stationary or non-stationary - is explicitly given by the factor 
exp(— Fti/2) in the amplitude ( p.2|) according to the Weisskopf-Wigner approximation. As 
shown in Section III the correction factor gi2 depends on the quantity 

p = ——4 , (4.2) 

and is only sizeably different from one if p > 1. Thus a convenient formulation of SFC is 
p < 1. Note that by defining a spread in velocity of the source wave packet by Avs = cts/^ij. 
the last formulation of SFC can be rewritten in analogy to Eq. (^.1|) as |M] 
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(4.3) 



As we have seen in Section III with as ~ 0.01 MeV one gets p ~ 0.005 and SFC is well 
fulfilled for the LSND and KARMEN experiments, however, the margin for violation of SFC 
is only two orders of magnitude, in contrast to ACC where the margin is at least ten orders 
of magnitude. 

Taking this observation as a starting point, we want to see if there is some chance to 
increase p to 1 or more in order to have an observable suppression effect. Looking at Fig. 2 
one recognizes that in the case of erg = 10 MeV the lower line of the curve has undergone a 
considerable shift downwards compared to the case p = 0. Note that this shift happens in 
the sensitive region of the LSND and KARMEN experiments. One might therefore ask the 
question if it is possible to increase as in order to achieve p ^ 1. If the source wave function 
is such that ascyxs ~ 1? which is correct for Gaussian wave functions, one would have to 
localize the neutrino source particle very well in coordinate space, namely to axs ^ 10~^^ 
m for as ^ 0.1 MeV. Taking a look at Eq. ( |4.2| ), one observes that p can also be increased 
by decreasing the neutrino energy. However, in this case one lowers the ratio of true over 
background events. Finally, a small product of mass times decay width of the source particle 
also enhances p but, unfortunately, the muon has already the smallest such product among 
the particles which can copiously be produced.0 

The quantum- mechanical ACC and SFC in the forms ( [4. 1| ) and ( [1.3| ), respectively, have 
analogous "classical macroscopic" conditions stemming from the inescapable averaging over 
some regions of the target and the detector, respectively, when an experiment is performed. 
Replacing in Eq. ( [4.1|) the widths axs,xD by the typical sizes Rs,d, respectively, of the regions 
in the target where the muon is stopped and in the detector where the neutrino detection 
process is localized we arrive at conditions which are obtained by the incoherent averaging 
of the oscillation probabilities over the variations of L due to Rs,d and the requirement that 
neutrino oscillations should not be washed out by this averaging process. Since Rs,d are 
macroscopic quantities, if these classical conditions are fulfilled, then clearly also ACC holds 
because axs,xD Rs,d- The classical analogue to SFC ( ^I3| ) says that during its lifetime 
the neutrino source particle should move a distance much less than the oscillation length in 
order not to wash out neutrino oscillations. This classical condition is not obviously linked 
to SFC. 



^The corresponding product for the neutron is smaller, however, to observe this effect in neutri- 
nos from neutron decay one would have to isolate the neutrons from the external interactions - 
something that seems to be hard to achieve in experiments. 
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FIGURES 
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FIG. 2. Contour plot for the transition probability P = J (f){E^)Pp^^p^ {E„)dE^ at P = 0.001 in 
the Am? ~ sin^ 26 plane where Pp^^^p^ is defined in Eq. (|3.4D and (p denotes the energy distribution 
of the flux (see, e.g., Ref. [16]). We integrated over the neutrino energy Ei, from 12 up to 53 
MeV and used L = 17.7 m for the source ~ detector distance L. The solid, dashed and dotted line 
correspond to as = {g = 1), 1 and 10 MeV, respectively. 
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